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ABSTRACT 

The kinematics of the extra-planar neutral and ionised gas in disc galaxies shows 
a systematic decline of the rotational velocity with height from the plane (vertical 
gradient). This feature is not expected for a barotropic gas, whilst it is well reproduced 
by baroclinic fluid homogeneous models. The problem with the latter is that they 
require gas temperatures (above 10^ K) much higher than the temperatures of the 
cold and warm components of the extra-planar gas layer. In this paper, we attempt 
to overcome this problem by describing the extra-planar gas as a system of gas clouds 
obeying the Jeans equations. In particular, we consider models having the observed 
extra-planar gas distribution and gravitational potential of the disc galaxy NGC891: 
for each model we construct pseudo-data cubes and we compare them with the HI 
data cube of NGC 891. In all cases the rotational velocity gradients are in qualitative 
agreement with the observations, but the synthetic and the observed data cubes of 
NGC 891 show systematic differences that cannot be accommodated by any of the 
explored models. We conclude that the extra-planar gas in disc galaxies cannot be 
satisfactorily described by a stationary Jeans-like system of gas clouds. 

Key words: galaxies: kinematics and dynamics - galaxies: haloes - galaxies: individ- 
ual: NGC 891 - galaxies: ISM - galaxies: structure - ISM: kinematics and dynamics 



1 INTRODUCTION 

Observations of several spiral galaxies at various wave- 
lengths have revealed massive gaseous haloes surrounding 
the galactic discs. This extra-planar gas is multiphase: it is 
detected in Hi (e.g. Swaters, Sancisi & van der Hulst 1997), 
in Ha (Rand 2000; Rossa et al. 2004) and in X-ray obser- 
vations (Wang et al. 2001; Strickland et al. 2004). Thanks 
to high-sensitivity Hi observations of some edge-on galax- 
ies, like NGC 891 (Oosterloo, Fraternali & Sancisi 2007), 
we can trace these haloes up to 10 - 20 kpc from the plane 
and we can study their kinematics in detail. A remarkable 
feature of the extra-planar gas is its regularly decreasing ro- 
tational velocity at increasing distances from the galactic 
plane (vertical gradient). Fraternali et al. (2005) found that 
in NGC 891 the vertical gradient is ^ — 15km s~^kpc~^ (see 
also Heald, Rand & Benjamin et al. 2007). For other galax- 
ies (e.g. NGC 2403) gradients of the same order have been 
estimated (see Fraternali 2008). 

The dynamical state of the extra-planar gas can provide 
useful insights to understand its origin. Two main frame- 
works have been proposed: the galactic fountain (Shapiro & 
Field 1976; Bregman 1980) and the cosmological accretion 
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(Oort 1970; Binney 2005). In the galactic fountain model, 
partially ionized gas is ejected from the disc by supernova 
explosions and stellar winds, it travels through the halo 
and eventually falls back to the disc. Due to difficulties in 
running high-resolution hydrodynamical simulations of the 
whole galactic disc (see e.g. Melioli et al. 2008), the gas 
in the galactic fountain model is often treated ballistically, 
i.e. as a collection of non-interacting clouds subject to the 
gravitational field of the galaxy (Collins, Benjamin & Rand 
2002). The ballistic galactic fountain is able to reproduce the 
gas distribution observed in spiral galaxies, but in general 
the predicted vertical gradient of is significantly lower 
than observed (Fraternali & Binney 2006). In the accretion 
models, the extra-planar gas is the result of the infall of 
cool intergalactic gas. The two scenarios are not mutually 
exclusive: in fact it has been argued that the contribution 
of accretion is limited to only 10 — 20% of the total mass of 
extra-planar gas (Fraternali & Binney 2008). 

In a complementary approach, some authors (Benjamin 
2002; Barnabe et al. 2006), have focused on the construction 
of models for the extra-planar gas under the assumption of 
stationarity. In this approach the problem of the origin of 
the extra-planar gas is not addressed, and the main effort is 
to clarify the global dynamical state of the gas. Usually, in 
the stationary approach, the extra-planar gas is described as 
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an homogeneous fluid in permanent rotation. In the simplest 
of fluid homogeneous models, vertical and radial motions of 
the gas are neglected and therefore the vertical gravitational 
field of the galaxy is balanced only by the pressure gradi- 
ent of the gas. In other studies, turbulence has also been 
added as an extra- pressure term (e.g., Koyama & Ostriker 
2008). According to the Poincare-Wavre theorem (Lebovitz 
1967; Tassoul 1978), a vertical gradient in the rotational ve- 
locity can be present only if the pressure of the medium 
is not stratified on the density, therefore in fluid homoge- 
neous models the gas distribution is necessarily baroclinic 
(Barnabe et al. 2006; see also Waxman 1978). Remarkably, 
in their study of NGC 891, Barnabe et al. (2006) found that 
if the baroclinic configuration is built from general physical 
arguments, the model rotational velocity well reproduces the 
observed rotation curves at different heights over the galactic 
disc and for a large radial range. However, the temperature 
predicted for the system is > 10^ K, well above that of the 
neutral gas. 

Here we investigate if and how this "temperature prob- 
lem" can be solved. We make use of the fact that stationary 
fluid equations for a gaseous system in permanent rotation 
are, from a formal point of view, identical to the stationary 
Jeans equations for an axisymmetric system with isotropic 
velocity dispersion. Thus, also in this case, a decrease of 
rotational velocities with increasing distances from the mid- 
plane is obtained when the density field is built using the 
approach described by Barnabe et al. (2006). In this paper 
we analyse this alternative interpretation of the baroclinic 
fluid homogeneous models. In practice we consider a "gas" 
of cold HI clouds described by the stationary Jeans equa- 
tions, where the thermal pressure, needed for the vertical 
balance of the galaxy gravitational field, is replaced by the 
velocity dispersion of the clouds. 

In the present approach, some preliminary considera- 
tion on the physical state of the clouds is important. In fact, 
Hi halo clouds are often assumed to be embedded in, and 
in pressure equilibrium with, a hot medium (corona) that is 
about a thousand times less dense and provides the pressure 
required to confine them (Spitzer 1956; Wolfire et al. 1995). 
In the Milky Way, there is plenty of evidence for the exis- 
tence of this medium, such as emission from highly ionised 
metals (e.g. Sembach et al. 2003) and the head-tail morphol- 
ogy of individual High Velocity Clouds (HVCs, e.g. Briins 
et al. 2000). Due to the high density contrast, the pressure- 
confined clouds cannot be supported by the pressure of the 
external medium and must be moving more or less ballisti- 
cally in a way akin to water drops in the air (Bregman 1980). 
This assumption is commonly employed in galactic fountain 
models (e.g. Collins et al. 2002). On the time-scale taken 
for a cloud to move about 1000 times its length, the trajec- 
tory may deviate significantly from the ballistic one due to 
the interactions with the external medium. This time-scale 
depends on the cloud mass but it is longer than a dynam- 
ical time for the larger clouds (Fraternali & Binney 2008). 
The typical mass of an halo cloud is difficult to estimate. 
In the Milky Way the new determinations of distances for 
the HVCs give masses from a few xlO^ M© to 10^ M© (e.g. 
Wakker et al. 2008), for external galaxies the mass resolu- 
tion is often above 10^ Mq but in M31 several clouds with 
masses down to ~ 1 x 10^ Mq have been observed (Thilker 
et al. 2004). We estimate that for a cloud of 10^ M©, the 



drag time scale is more than 10 times the dynamical time 
and it is therefore fair to treat the system as purely dy- 
namical. Moreover, if the HI halo of a galaxy like NGC 891 
is made up of clouds with this mass and radii of 100 pc, 
the collision time between clouds turns out to be about 5 
times the dynamical time and the system can be considered 
collision-less. 

The substitution of a fluid system with a gas of clouds is 
also not trivial from the point of view of the comparison with 
the observations, as several delicate issues arise when con- 
sidering Hi observations (see Sect. 13. 3|) . To overcome these 
problems we follow Fraternali & Binney (2006, 2008) and 
construct, as an output for each of the models investigated, 
a pseudo-data cube with the same resolution and total flux as 
the data cube of the Hi observations. This procedure assures 
a full control of projection and resolution effects. Moreover, 
the comparison with the raw data removes all the intermedi- 
ate stages of data analysis and the associated uncertainties. 
For completeness, we also investigate some phenomenologi- 
cal anisotropic models. 

The paper is organized as follows: in Section 2 we illus- 
trate the Jeans-based interpretation of the baroclinic solu- 
tions, introducing the anisotropic case which has no anal- 
ogous fluid counterpart, and we discuss the conditions re- 
quired for these solutions to have a negative vertical veloc- 
ity gradient. In Section 3 we present the method adopted to 
construct isotropic and anisotropic models for the edge-on 
galaxy NGC 891, whilst in Section 4 we compare their pre- 
dictions with the Hi observations. Section 5 is devoted to 
the discussion of the results, and Section 6 concludes. 



2 JEANS-BASED DESCRIPTION OF FLUID 
BAROCLINIC SOLUTIONS 

2.1 Isotropic Jeans solutions 

Consider an axisymmetric density distribution of gas clouds 
moving under the influence of an axisymmetric galactic grav- 
itational potential z), so that all the physical quanti- 
ties are independent of the azimuthal angle cp of cylindrical 
coordinates. The galactic gravitational potential $ is the 
sum of the dark-matter halo, the bulge and the stellar and 
gaseous disc potentials. We neglect the contribution to the 
gravitational field of the extra-planar gas, which represents 
a very small fraction of the total mass (for instance, less 
than 1% in NGC 891). If the velocity dispersion tensor of 
the cloud distribution is isotropic the associated stationary 
Jeans equations are 

p dz 

_ (1) 

1 dpal ^ yI - al 

p dR dR R ' 

where p is the cloud density distribution, an = cr^ = a is the 



dz 



one- dimensional velocity dispersion, and is the velocity 
of each cloud around the axis (e.g., Binney & Tremaine 
2008). The phase-space average is indicated by a bar over 
the quantity of interest so = V(^, and we are assuming 
that ur = Uz = 0, i.e. no net motion in the radial and the 
vertical directions is considered. Under these assumptions 
the streaming velocity of the clouds is given by 



Stationary models for the extra-planar gas in disc galaxies 3 



2 ~2" 2 



(2) 



and eqs. ([T]) are formally identical to the equations describ- 
ing a stationary fluid in permanent rotation, where the ther- 
modynamic pressure of the fluid is replaced by the quantity 
pa'. 

To solve eqs. {[1]) we fix the galactic gravitational po- 
tential ^{R,z) and we assume a cloud density distribution 
p{R,z) consistent with the observations and vanishing at in- 
finity (see Section[3]). With this procedure, the obtained con- 
figuration is in general baroclinic. Therefore, from the first 
of eqs. ([!]) with boundary condition pa^ for z ^ oo, 
one obtains 



per, 



oz 



(3) 



and is calculated from the second of the eqs. ([T]) and 
eq. 0. Note that the rotational velocity field can be 



obtained in general by the "commutator-like" relation 



R 



[or dz 



dp d^' 
d^dR. 



dz , 



(4) 



(e.g., see Barnabe et al. 2006; and references therein). 

The expression above reveals that the existence of phys- 
ically acceptable solutions with ^ everywhere is not 
guaranteed for independent choices of p and In Barnabe 
et al. (2006) general rules for constructing physically ac- 
ceptable baroclinic solutions are illustrated. For example, 
in the usual case in which d^/dR ^ and d^/dz ^ 0, 
a sufficient condition to ensure the positivity of is that 
dp/dz ^ and dp/dR ^ 0. The same rules apply also here 
suggesting the use of centrally depressed density distribu- 
tions. This implies that in the central regions the density 
is radially increasing moving outwards, which is consistent 
with the radial Hi distribution observed in spiral galaxies 
such as NGC 891 (Oosterloo et al. 2007). 



2.2 Anisotropic Jeans solutions 

In Section [3] we will show that Jeans isotropic models, the 
natural kinetic counterpart of homogeneous fluid models, 
fail to reproduce some features of the extra-planar gas in 
NGC 891. For this reason we also explore a phenomenolog- 
ical model including velocity dispersion anisotropy. Given 
the gravitational potential $ and the cloud density distri- 
bution p, the solution of the stationary axisymmetric Jeans 
equations with an anisotropic velocity dispersion tensor is 
unique if one assigns the orientation and the shape of the 
intersection of the velocity dispersion ellipsoid everywhere in 
the meridional plane. For our problem, this is done by main- 
taining ur = Uz = and assuming a suitable parametriza- 
tion for the anisotropy, which specifies the ratio grIgz - The 
simplest choice is that of a constant anisotropy which, fol- 
lowing Cappellari (2008), we express as a\ hal with 6^0 
0. Thus, eqs. ([T]) become 



^ Note that this anisotropy cannot be realized (except for h = 
1) under the common assumption of a two-integral distribution 
function depending on energy and angular momentum along the 
2;— axis, as in this case (Tr = Gz- 
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Thus erf is given by eq. JS]) as in the isotropic case, 
whilst it can be shown that 
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Note that equation above can be rewritten in a more com- 
pact form as 



(7) 



where Uiso is the rotational velocity field of the correspond- 
ing isotropic model (defined as the model with the same 
potential and the same density distribution, but with k = 
6=1) and 



R 



OR 



(8) 



is the circular velocity at position {R,z). The streaming ve- 
locity field is obtained by splitting vj into the streaming 
motion and the azimuthal velocity dispersion crj. Fol- 
lowing the widely used Satoh (1980) /c- decomposition we 
have 



ul, = /c^(vj - baz) 



so that 



cr^ = ba^ + (1 - (vi 



(9) 



(10) 



with ^ /c ^ 1; in general k can be function of the position, 
i.e. k = k{R, z) (e.g., see Ciotti & Pellegrini 1996). For /c = 
no ordered motions are present, whilst a J = cr^ in the case 
A; = 1; the isotropic case is recovered for k — b — 1. From eq. 
([9|) it follows that the Satoh decomposition procedure can 
be applied only when v^ — 6crf ^ 0, so /c and b are not fully 
independent. Note that if the corresponding isotropic model 
leads to a physically acceptable solution, then a sufficient 
condition to ensure the positivity oi is b ^ 1\ \i u(^^ ^ 0, 
the l.h.s. of eq. (0) can be negative only if 6 > 1, so the use 
of the Satoh decomposition is consistent with the assump- 
tion gr ^ Gz. This requirement is in agreement with the 
observations that show large vertical motions of the order 
of 50 — lOOkms"^ and velocity dispersion along R and of 
the extra-planar gas of approximately 20kms~^ (Boomsma 
et al. 2008; Oosterloo et al. 2007). Summarizing, this simple 
parametrization of the velocity dispersion anisotropy is con- 
sistent with physically acceptable solutions provided that: 
(i) 6 ^ 1 and (ii) the corresponding isotropic model has 
ut,^ everywhere. 



2.3 Conditions to have negative vertical gradient 

An important property of the extra-planar gas is the de- 
crease of its rotational velocity with the distance from the 
galactic plane. To study the kinematics of the extra-planar 
gas in the models it is helpful to have general expressions 
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(such as eq. [2] of Waxman 1978) relating the vertical gra- 
dient of the rotational velocity to other physical quantities 
in stationary configurations. 

For example, in the isotropic case we differentiate eq. 
Q with respect to z and obtain 



d\np 



dz dz ^isoj d\nR dz 

Thus, the necessary and sufficient condition to have a de- 
crease in the rotational velocity in the isotropic model is the 
negativity of the right hand side of eq. In particular, 

consider a typical potential such that /dz ^ and a cen- 
trally depressed density distribution vanishing for ^ oo 
with dp/dz ^ 0. For each value of z there is a radius i?max 
such that p increases with R ioi R < Rmax and then de- 
creases with R for R ^ i?max- From eq. it follows that 
for R ^ i?max a necessary condition for a solution to ex- 
hibit a negative vertical gradient in the rotational velocity 

IS IXiso ^ Uc . 

For systems with anisotropic velocity dispersion tensor 
of the class here considered, there is not a simple expression 
for the vertical gradient of the rotational velocity, unless we 
assume that k is independent of z. In this case, a useful 
expression for the vertical gradient of the rotational velocity 
can be obtained by substituting eq. (|9]) into eq. (0 and by 
differentiating the result with respect to z. This leads to 



d\np 



(11) 



dz 



dz dz 



2 

iso 



(12) 



where i^iso is the rotational velocity of the corresponding 
isotropic model. It follows that in the usual cases, in which 
du^/dz ^ and b < 1 (see Section [2.2p . a negative verti- 
cal velocity gradient in the corresponding isotropic model 
is a sufficient condition to have du^jdz ^ also in the 
anisotropic models of the family here considered. 



3 APPLICATION TO NGC 891 

The general rules of Section 2 are now applied to the specific 
case of the extra-planar gas of NGC 891, to explore whether 
the proposed models reproduce the observed decrease of the 
rotational velocity for increasing z. We note that the result- 
ing models can directly describe the extra-planar HI gas seen 
in the observations because the vertical "pressure" support, 
needed against the galaxy gravitational field, is now pro- 
vided by the velocity dispersion of the clouds, rather than 
by the thermodynamic temperature as in the homogeneous 
baroclinic models of Barnabe et al. (2006). The present in- 
vestigation differs from that of Barnabe et al. (2006) also be- 
cause we adopt an observationally justified density distribu- 
tion for the extra-planar gas and we perform the comparison 
of the models with the observational data by constructing a 
pseudo-data cube for each model with the same resolution 
and total flux as the observations. 



3.1 The galaxy model 

For the determination of the gravitational potential of NGC 
891 we consider the maximum- light mass model of Fraternali 
& Binney (2006). This mass model consists of four distinct 
components: two exponential discs (stellar and gaseous), a 
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Figure 1. Total HI emission map of NGC 891 as obtained by 
the observations (red contours and colour scale) and by the cloud 
density distribution of Section [221 (black contours). The map has 
a spatial resolution of 25^^ (~ 1.2 kpc) with contour levels, starting 
from the outermost, of 0.15, 0.3, 0.6, 1.2, 2.4, 4.8, 9.6 xlO^o cm-^. 
The shaded area delimits the emission expected from the disc. 



bulge and a dark-matter halo. The volume density of the 
two disc components is given by 



p(7?,0) = poe-^/^<isech2 (^) , 



(13) 



where po is the central density, and the dimensional constant 
Rd and hd are, respectively, the scale-length and the scale- 
height of the disc. The bulge and the dark halo components 
are modelled using the double power-law spheroidal density 
profile (Dehnen & Binney 1998) 



p(m) 



where 



(mjaY (1 + m/a)^- 



m — R ^ 



(14) 



(15) 



e is the eccentricity of the spheroid, a and /3 are the inner 
and outer slopes, and a is the scale-length. The adopted 
values of the parameters to reproduce the rotation curve in 
the plane of NGC 891 are reported in Table 1 of Fraternali 
& Binney (2006), where the derivation of this potential and 
other mass models compatible with this rotation curve are 
also discussed. 



3.2 The cloud distribution 

A proper Jeans-based analysis of baroclinic fluid solutions 
requires the choice of a cloud density distribution similar to 
that observed. For the cloud density distribution we adopt 
the function 



p{R, 



+ — 

Rq- 



-R/ Rg 



sech -— 



(16) 



which is similar to that adopted by Oosterloo et al. (2007). 
In the equation above po is the central density, 7 ^ 0, i?g is 
the scale-length and the scale-height hg depends upon R as 
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h^{R) = ho+ , S^o. (17) 

The parameters of the distribution, excluded the central 
density po, are chosen is such a way that the function re- 
produces the Hi total map shown in Fig. [T]and their values 
are: Rg = 1.61 kpc, 7 = 4.87, ho = 2.3 kpc, 6 = 2.25 and 
Hr = 18.03 kpc. The central density po = 6.63 x 10~^cm~^ 
is such that the total mass of the cloud system is that de- 
rived by the total Hi emission; the corresponding halo mass 
(i.e. the mass computed for \z\ ^ 1 kpc) of 1.13 x lO^MO is 
consistent with the value given in Oosterloo et al. (2007). 

We now verify that the density distribution guarantees 
the existence of physically acceptable solutions (i.e. 
everywhere). It can be easily shown that dp/dz ^ for the 
density distribution {[16}, whilst the radial density profile 
presents, for each a maximum at some value Rma^{z): the 
positivity of is then guaranteed in the region R ^ i?max, 
where dp/dR > 0, whilst it has to be investigated in the ex- 
ternal region, where dp/dR < (see Sect ion [24]), by numeri- 
cal inspection. We found that, with the adopted parameters, 
the density distribution ([T6)) produces physically acceptable 
solutions for all the models investigated. For what concerns 
the vertical velocity gradient, the considerations made in 
Section 2.3 are valid: in particular, we expect that the ro- 
tational velocities of the models are lower than the circular 
velocity {[8]) for ^ i?max, if the vertical velocity gradient is 
negative. 

In our models, we are only aiming to reproduce the 
extra-planar gas and we therefore neglect the region where 
emission from the disc is expected. Such a region is cov- 
ered by shaded areas in all the figures where the comparison 
with the data is shown (Figs. [1] [3] [5] and [6}. The size of 
this region {\z\ < 1.64 kpc), is calculated as ±3cr x (spa- 
tial resolution of the data) assuming that the galactic disc 
is perfectly^ co-planar. In fact the resolution of the observa- 
tions (HPBW= 0.64 kpc) is much larger than the expected 
intrinsic scale height of the Hi disc (100 — 200 pc), i.e. the 
disc is not resolved. The exclusion of the disc is necessary to 
avoid the complication of a region where the gas is almost 
certainly not described by the considered Jeans equations. 
The implicit assumption is that, in a real galaxy, at the 
lower halo boundaries the velocity dispersion, induced by 
some (unspecified) physical process in the disc, is the same 
as the one provided by our density distribution (see Figs. [2] 
and [3]). This requirement is necessary for the system to be 
stationary. 



3.3 Construction of the model 

Having fixed the galaxy gravitational field and the cloud 
density distribution, for each model we first compute az 
by integrating numerically eq. ([3]) with a finite- difference 
scheme. In the isotropic case we then compute the rotational 
field from eq. ([4]) using the same numerical scheme. In the 
anisotropic cases we can integrate eq. ([6]) or just use ([T]). 
The rotation velocity and the azimuthal component of the 
velocity dispersion are then evaluated on the numerical grid 
by means of eqs. ([9]) and (|10|) . We verified that eqs. ([6]) and 
give the same results in all explored anisotropic models. 

As an output for the models a pseudo-data cube, with 
the same spectral and spatial resolutions and total flux as 



the data cube of Hi observations, is generated as follows. 
In each grid point the cloud density, the rotational velocity 
and the velocity dispersion are projected numerically along 
the line-of- sight for the inclination of the galaxy, which is 
assumed to be seen exactly edge-on. From the projected ve- 
locity we construct a first data cube and, to account for the 
effect of the velocity dispersion, we convolve the projected 
density in each cell of this cube with a normalised Gaus- 
sian profile with a dispersion equal to the projected veloc- 
ity dispersion. Finally, for a more realistic comparison, the 
cube is smoothed to the same spatial resolution as the data 
(FPBW = 28" ^1.3 kpc). We remind that the construc- 
tion of a pseudo-data cube is crucial because it bypasses all 
the intermediate stages of data analysis, and the assump- 
tions made in these stages. A critical case is the comparison 
of the data rotation curves with the rotational velocities pre- 
dicted by the models. The rotation curves of the extra-planar 
gas (and of edge-on galaxies in general), are usually derived 
by fixing a certain amount of line-of-sight velocity dispersion 
and using the so-called "envelope tracing" method (see San- 
cisi & Allen 1979; Fraternah et al. 2005). Clearly, if a model 
has velocity dispersions inconsistent with that adopted to 
derive the rotation curves from the data (as it happens here, 
see Section [4} , a direct comparison between the two cannot 
be done (see Section [5}. 



4 RESULTS 

We now present the most important features of some of 
the considered stationary models for the extra-planar gas 
of NGC 891, and we compare them with the observations. 

4.1 Isotropic model 

We start by describing the main intrinsic features of the 
isotropic model, i.e. the analog of a stationary baroclinic 
fluid model where the thermal pressure of the gas is replaced 
by a globally isotropic velocity dispersion tensor. It is then 
expected that the behaviour of the isotropic model is similar 
to the baroclinic solutions of Barnabe et al. (2006). However, 
there is an important difference between these two models in 
the choice of the density distribution. In our model the den- 
sity is chosen to match the observed Hi emission (see Fig. 
[U, whilst the density distribution of Barnabe et al. (2006) 
was just chosen to satisfy the general rules that guaran- 
tee a physically acceptable rotational velocity field. In other 
words, our density distribution is fixed by the Hi observa- 
tions, and not by the requirement ^ 0. Using this cloud 
distribution the resulting az = ctr = has values ranging 
from 160 (in the centre) to 35kms~^ (at R ^ 20 kpc) with 
mean values between 50 and lOOkms"^ (Fig. [2] bottom left 
panel) . 

In the left panels of Fig.[2]we present the rotation curves 
(top) and the radial velocity dispersion profiles (bottom) at 
distances from the plane 2; = 1.7 kpc (lower halo boundary) 
and z = 5.1 kpc; the dashed lines show the circular velocities 
Uc. As expected from the considerations of Section 2.3, in 
the external regions of the galaxy the rotational velocity is 
below the circular velocity because of the negative vertical 
gradient in u^. The radial profile of the velocity dispersion 
shows a monotonic decrease with R but the values of the 
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Figure 2. Intrinsic rotational velocities and velocity dispersion for system of clouds that can be described by a classical two-integral 
phase-space distribution function, at 2: = 1.7 kpc and z = 5.1 kpc. From left to right: isotropic model, k = 0.9 model and maximally 
rotating model (a^p = 0). Velocity panels: the dotted lines represent the circular velocity Uc whilst the solid lines are the model rotational 
velocities u^p. Dispersion panels: the dashed lines indicate az and the solid lines a^p. In the isotropic case all the components of the 
dispersion are the same, whilst in the other two models az = ctr. 



dispersion are everywhere greater than the mean observed 
value of 20kms~^. This latter is also the value adopted 
to derive the rotation curves above the plane from the Hi 
observations of NGC891 (Fraternali et al. 2005, here shown 
in Fig. [T]). As a consequence, a direct comparison between 
these curves and the circular velocities of our models is not 
possible. We stress again that these problems are overcome 
by comparing directly the models with the raw data cubes, 
where all the kinematical and density effects are fully taken 
into account. 

Figure [3] (first column from the left) shows the pseudo- 
data cube constructed from this Jeans model, while each 
panel in the rightmost column represents the observed emis- 
sion of Hi gas in a particular velocity channel. The numbers 
on the top right corners of the data column represent the 
deviations of the line-of-sight velocities from the systemic 
one (528kms~^), whilst the shaded area on each map de- 
limits the emission expected from the disc (zb3cr spatial res- 
olution), which we are not attempting to reproduce with 
this model. We can see that the observed Hi emission at 
higher rotational velocity (corresponding to the channel at 
'^^los = — 236.0kms~^) is fully due to the mid-plane, i.e. the 
halo is not visible. This is not surprising because the halo 
gas rotates more slowly than the gas in the disc and in fact, 
as we move towards the bottom panels, the halo emission 
becomes important. Due to its high velocity dispersion, in 



the isotropic Jeans model some emission in the halo region is 
present also in the upper two velocity channels; only in the 
middle channel the agreement is good whilst in the bottom 
channels the emission is too horizontally elongated. 

Summarizing, it is possible to construct an isotropic 
Jeans model with a realistic halo density distribution, in 
which a negative vertical gradient in the rotational veloc- 
ity is present. However, the predicted velocity dispersions 
along the line-of-sight are much higher than those observed, 
so that a number of features of the data cube are not repro- 
duced. In order to solve these problems, in the next Section 
we consider different models with anisotropic velocity dis- 
persion tensors. 



4.2 Anisotropic models 

The assumption of anisotropy in the velocity dispersion can 
overcome, at least in principle, the problems originated by 
the high values of a obtained in the isotropic models. In 
Section 14.2.1 1 we focus on the classical Satoh decomposition 
in a two- integral distribution function (i.e., gr — Gz every- 
where) , and in Section 14.2.21 we discuss the case of a three- 
integral phase-space distribution function for the cloud sys- 
tem. In the following, we indicate with aios the component of 
the velocity dispersion along the line-of-sight on each place 
in the galaxy. 
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Figure 3. Comparison of representative velocity channels of the HI data cube of NGC 891 and the predictions of 3 stationary models: 
the isotropic (first column), the k = 0.9 and 6 = 1 (second column), and the b = 1 maximally rotating (third column) models. The 
shaded area in each channel delimits the emission expected from the disc. The contour levels, starting from the outer, are 0.25, 0.5, 1, 
2, 5, 10, 20, 50 mJy/beam; the beam size is 28''= 1.3 kpc. 



4-2.1 Two-integral anisotropic systems 

We examine here cloud systems described by a two-integral 
phase-space distribution function, for which az = o-r. The 
relevant equations are obtained by taking 6 = 1 in eqs. (|5|), 
and adopting a Satoh decomposition in the azimuthal di- 



rection. When /c = crj is maximal, while when /c = 1 
the system is isotropic. Allowing for k = k{R,z), more ro- 
tationally supported models can be constructed, up to the 
maximally rotating one where cr^ = (Ciotti & Pellegrini 
1996). First of all, we investigate a model with constant 
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Figure 4. Intrinsic rotational velocities and velocity dispersion for the anisotropic model 6 = 0.1 and k = 1 (left column), the model 
b = 0.3 and k = 1 (central column), and the model 6 = 0.1 and k = 0.9 (right column) at z = 1.7 kpc and z = 5.1 kpc. Rotational 
velocity panels: the dotted lines represent the circular velocity Uc whilst the solid lines are the rotational velocities u,p. Dispersion panels: 
the dashed lines indicate cr^, the solid lines a^p and the dashed-dotted lines aR] in the first two cases a^p = ctr. 



k and we assume k = 0.9; for lower values the disagree- 
ment with the data increases. We expect that the model is 
not better than the isotropic (/c = 1) case in reproducing 
the observations, because aios is increased by the adopted 
anisotropy; in particular, as can be inferred from eq. (|10|) , 
(Jip ^ cTz = an. The central column of Fig. [2] (bottom panel) 
clearly illustrates this fact, showing a difference among 
and the other two components of approximately 30kms~^, 
on average. For what concerns the properties of the rotation 
curves, the considerations carried out in the isotropic case 
are still valid, because from eqs. (0) and (|9]) it follows that 
= kuiso so that in the external regions of the galaxy 
< Uc. The pseudo-data cube of the model presented in 
Fig. [3] (second column) shows that in the channels close to 
the systemic velocity (bottom) the emission of this model is 
even more elongated than in the isotropic case, in complete 
disagreement with the data. 

In two-integral systems, the line-of-sight velocity dis- 
persion can be reduced only if k > 1. As a limit case, 
we present the results for the maximally rotating model in 
which k = kmax{R,z) ^ 1 such that the rotational velocity 
has its maximum value everywhere and = 0. The cor- 
responding intrinsic rotation curves are presented in Fig. [2] 
(top right panel) and their peculiarity is the steep rise in 
the central regions. The reduced value of aios improves the 
agreement with the data, as it can be seen from the pseudo- 
data cube of the model in Fig. [3l especially in the channels 



near the systemic velocity. This agreement is not so good 
in the channel at v\os = — 236.0kms~^ and in general in 
all the channels between —236.0 km s~^ and —161.8 km s~^ 
(not shown here) where, due to the steep rise in the rota- 
tional velocity at small R, there is emission in the central 
halo regions not seen in the data. 



4-2.2 Three-integral anisotropic systems 

We now turn the attention to the models with b ^ 1 (there- 
fore supported by a three-integral phase-space distribution 
function) and we first discuss the case an — (i.e. k — 1). 
The best fit to the data is obtained by an anisotropic model 
with b ^ 0.1, which limits the line-of-sight velocity disper- 
sion to about 20kms~^ for a large interval of R (Fig. 21 
bottom left panel). The rotation curves for this model (top 
left panel) are very similar to the circular velocity curves, 
in accordance with eq. (0 which states that, in the limit- 
ing case 6 = 0, the curves coincide. The pseudo-data cube 
of the model is presented in Fig. [5] (leftmost column). In 
the high velocity channels (top panels) the halo emission is 
still present but, at variance with the isotropic case, this is 
due to the fact that the rotational velocities of the model 
are slightly higher than those observed. In the low velocity 
channels (bottom panels) the halo emission is narrower than 
the data along the horizontal axis, due to a rapid decrease of 
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Figure 5. Comparison of five representative velocity channels of the HI data cube of NGC 891 and the predictions of 3 three-integral 
anisotropic models: 6 = 0.1 and k = 1 (first column), b = 0.3 and k = 1 (second column), 6 = 0.1 and k = 0.9 (third column). The 
shaded area in each channel delimits the emission expected from the disc and the contour levels are the same as in Fig. [3] 



the velocity dispersion in the outer parts. We discuss these 
discrepancies more in detail in Section [5] 

An important issue regarding the anisotropic models is 
the effect that the choice of the parameters b and k may 
have on the properties of the solutions. We analyse this as- 
pect by fixing the value of one of the two parameters and by 



varying the other. The results of this comparison, for rep- 
resentative values of b and k, are shown in Fig. |4l where 
the model with 6 = 0.1 and /c = 1 (which, hereafter, we 
will also refer to as the "high-cr^ model") has been taken 
as a reference case. In particular, increasing the value of b 
at fixed k has the net result of increasing the line-of-sight 
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velocity dispersion (i.e. gr and a,^) whilst the rotational ve- 
locity is decreased; the decrease in rotational velocity is also 
obtained when the value of k is decreased at fixed 6, but 
in this case only the value of is affected. In Fig. [5] we 
present the pseudo-data cubes of the models just described. 
The emission of the reference model (6 = 0.1, /c = 1) in all 
the velocity channels is less horizontally extended than that 
of the model with h = 0.3, due to the smaller values of aios- 
Although the 5 = 0.3 model seems to agree better with the 
data in the channel maps close to systemic velocity the top 
panels of Fig. [5] show substantial halo emission, due to the 
now increased line-of-sight velocity dispersion. The interpre- 
tation of the differences between the reference model and the 
model with /c / 1 is complicated by the fact that / gr 
so that the extension of the emission strongly depends on 
their relative contribution to the velocity dispersion along 
the line-of-sight. In particular, the importance of a,^ is max- 
imum in the high velocity channels (top panels) , and then it 
progressively diminishes in the channels near the systemic 
velocity in favour of gr. 

In short, three- integral models with b < 1 and k — 1 
perform better than the isotropic and the two-integral mod- 
els because cf\os can be reduced to values comparable to the 
the observed one. However, a number of features of the data 
cubes are still not reproduced revealing intrinsic problems 
with all these stationary models that we discuss more in 
detail in the next Sections. 

4.3 Position-velocity plots 

Among the models presented above we select the three most 
representative and compare them further with the HI data. 
Figure [6] shows two cuts of data cube of NGC 891 (rightmost 
column) taken perpendicularly to the plane of the galaxy at 
R = 7.5 kpc (first row) and R = 2.8 kpc (third row) from 
the centre, towards North-East. These plots are analogous 
to those shown in Fig. 16 of Oosterloo et al. (2007). The 
position- velocity diagrams are compared with the prediction 
of three models: the isotropic model, the maximum rotation 
model and the three- integral high-cr^^ model. Also the cor- 
responding residuals are shown (second and fourth rows) 
computed as the difference between the data and the mod- 
els; the shaded area in each posit ion- velocity plot delimits 
the emission expected from the disc. 

The characteristic triangular shape of the data dia- 
grams (rightmost column) is due to the negative vertical 
gradient in the rotational velocity of the halo and is the key 
feature that the models should reproduce. Clearly, this is 
not the case for the isotropic model as the high velocity dis- 
persion makes the emission very elongated, although a neg- 
ative vertical gradient in the rotational velocity is visible. 
In the 2.8 kpc cut (third row), the emission is so elongated 
that a substantial amount of gas is apparently "counter- 
rotating" (upper part of the plot). The maximally rotating 
model (second column) also does not reproduce the shape 
of the emission especially at = 2.8 kpc from the centre of 
the galaxy (third row). Also in this case the predicted line- 
of sight velocity dispersion is too high. The position- velocity 
diagrams for the high-cr^ model (third column) show a qual- 
itative agreement with the data as the triangular shape is 
reproduced. However the velocity dispersion in the upper 
halo > 5 kpc is clearly not matched and there is a down- 



ward shift of the emission especially visible in the 7.5 kpc 
cut. These features are also apparent in the residual maps 
and in particular the high-cr^ model, despite reproducing 
well the shape of the data, shows substantial residuals. In 
addition, the residuals also reveal the intrinsic asymmetric 
distribution of the extra-planar gas with respect to the mid- 
plane of NGC 891. 

Figure [6] gives a different view of the discrepancies al- 
ready noted in Figs. [3] and [5] The advantage of a cut per- 
pendicular to the plane of the galaxy is that it captures in 
one diagram the kinematics of the extra-planar gas at differ- 
ent heights and, as mentioned, its characteristic triangular 
shape is the evidence of the gradient in rotational velocity 
with height. The only model that reproduces at least quali- 
tatively this shape is the anisotropic 6 = 0.1, /c = 1 model. 



5 DISCUSSION 

In this paper we have explored an alternative interpretation 
of the baroclinic fluid homogeneous models for extra-planar 
gas, i.e. that of a "gas" of HI clouds described by the sta- 
tionary Jeans equations. Due to the formal equivalence of 
the equations describing the two systems, a similarity in 
the behaviour of the isotropic Jeans model and that of the 
baroclinic solutions of Barnabe et al. (2006) is expected. The 
present investigation is not just a re-interpretation of the re- 
sults of Barnabe et al. (2006): even in the isotropic case our 
models differ from theirs in the choice of galaxy mass and 
cloud density distributions. We find that the choice of the 
mass model has a little impact on the properties of the so- 
lutions, but the form of the density distribution is critical. 
Barnabe et al. (2006) adopted a very flat density distribu- 
tion, in order to reproduce the off-plane rotation curves. 
A flat density distribution requires less pressure support 
against the vertical gravitational field of the galaxy which 
results in a low temperature of the gas, typically 10^ K. 
For such a temperature the expected velocity dispersion in 
a Jeans isotropic model would have been of the order of 
30kms~^, close to the mean observed value (Oosterloo et 
al. 2007). Instead, in the present model the parameters of 
the density distribution are tuned to reproduce the observed 
HI distribution which has a greater scale-height than that 
used by Barnabe et al. (2006, see eq. HZ])- Thus, the typi- 
cal values of the velocity dispersion necessary to sustain the 
halo are of a factor 3 higher than the above. As we have 
seen, this high velocity dispersion has a large impact on the 
observable kinematics of the extra-planar gas and it is the 
main reason why the isotropic model fails to reproduce the 
data. The higher values of the velocity dispersion also in- 
fluence the rotational velocities: indeed the rotation curves 
of our Jeans model are about 20kms~^ lower than those 
obtained by Barnabe et al. (2006). 

The high values of the velocity dispersion of the Hi halo 
clouds may appear, at first sight, unphysical. However we 
notice that these velocities are unavoidable in any model at- 
tempting to reproduce gaseous halos with such large scale- 
heights (e.g. Fraternali & Binney 2006). The source of ki- 
netic energy for these clouds must come from supernovae, 
indeed hydrodynamical simulations of superbubble expan- 
sion in galactic discs predict these kinds of blowout veloci- 
ties (e.g. Mac Low & McCray 1988). Finally, we also notice 
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Figure 6. Comparison between two position-velocity cuts perpendicular to the plane of NCC 891 and some representative models. In 
the first row the cut is made at R = 7.5 kpc on side of the galaxy where the gas is moving towards the observer, and the corresponding 
residuals, are presented in the second row. In the third row the cut is on the same side of the galaxy but at J? = 2.8 kpc and the residuals 
are shown in the fourth row. The shaded area in each plot delimits the emission expected from the disc and the contour levels are the 
same as in Fig. [3] For the residuals, computed subtracting the model to the data, the black contours represents positive values whilst 
the white contours are the negative ones; also in this case the contour levels (absolute value) are the same as in Fig. [3] 



that velocities of the order of lOOkms"^ are in fact typical 
for halo clouds in the Milky Way (Wakker & van Woerden 
1997) and have been observed in external galaxies whenever 
the inclination of the disc along the line-of-sight is favorable 
(e.g. van der Hulst & Sancisi 1988; Boomsma et al. 2008). 

As already mentioned in Sect. 14.31 the most satisfac- 
tory model is the anisotropic 6 = 0.1, /c = 1 model (or 



high-cr^ model), which has been constructed in order to have 
ctr = < (Tz- These requirements are roughly in agreement 
with the prescriptions for a galactic fountain where the ma- 
terial is ejected mostly vertically from the disc and falls back 
also vertically. Fraternali & Binney (2006, 2008) have build 
galactic fountain models for NGC 891 and compared them 
with the data cube in the same way we did here. They found 
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that the halo gas distribution is well reproduced by fountain 
clouds that are kicked out of the disc with a Gaussian distri- 
bution of initial velocities and dispersion of about 80 kms~^, 
which is a value very similar to the average Gz predicted by 
our stationary models (Fig. |4]). Moreover, Fraternali & Bin- 
ney (2006) investigated the opening angles about the normal 
to the plane of the ejected clouds finding that it has to be 
lower than 15 degrees, in other words the clouds are ejected 
almost vertically. Therefore, the h — 0.1, /c = 1 model is 
quantitatively similar to the fountain models. What makes 
the difference between the two is the lack of non-circular (in- 
flow and outflow) ordered motions in our stationary model 
that, instead, in fountain models are always present. 

The most problematic discrepancy between the high- 
Gz model and the data is visible in the channel map close to 
the systemic velocity (see channel at — 5.2kms~^ in Fig. [5}. 
This map shows the gas that has a null velocity component 
along the line-of-sight (once the systemic velocity has been 
removed). Much of this gas is located in the outer parts (in 
R) of the galaxy both in the disc and in the halo. The width 
of the channel (along the horizontal axis) is mostly due to 
the velocity dispersion of this gas. The high- cr^ model fails 
to reproduce the data in two ways: 1) the channel map is 
too narrow, 2) the width (in the x-axis) of the emission 
decreases with height contrary to what seen in the data. 
These differences are caused by the fact that the velocity 
dispersion in the model decreases both with R and with 
whilst in the data the two behaviours seem reverse. The 
fact that the velocity dispersion increases with height has 
been clearly shown by Oosterloo et al. (2007) who built ar- 
tificial cubes with different parameters and compared them 
with the NGC 891 data cubes. These authors also kept the 
dispersion constant with R and reproduced much better the 
shape of the channel maps close to the systemic velocity (see 
their Fig. 14). The decrease of the velocity dispersion with 
R and 2; is a common problem of all our stationary models 
and may be related to the lack of (ordered) non-circular 
motions mentioned above. 

In the high-cr^ model the rotational velocities can be 
directly compared with those obtained from the data be- 
cause the line-of-sight velocity dispersion has values compa- 
rable with the observed one. In Fig. [7] we show the rotation 
curves of NGC 891 at distances from the plane of 3.9 kpc 
(top panel) and 5.2 kpc (from Fraternali et al. 2005). The 
rotational velocities predicted by the 6 = 0.1,/c = 1 model 
are shown as solid lines, the dashed line shows the rotation 
curve in the plane. In the derivation of these rotation curves 
from the data, it has been assumed a constant value of the 
velocity dispersiorEl (20kms"^). The curves have been de- 
rived using the envelope tracing method, which is such that 
the larger is the velocity dispersion one assumes the lower is 
the obtained rotational velocity. In our high- cr^ model the 
velocity dispersion is not constant along R and it is exactly 



^ Fraternali et al. (2005) assumed CTgas ~ Skms"-*^, for the disc 
and the halo. Oosterloo et al. (2007) found that the velocity dis- 
persion in the halo is actually larger and about 20kms~-'^, thus it 
is appropriate to shift the data points of Fraternali et al. (2005) 
down by llkms"-*^ (see Fraternali &: Binney 2008 for details). 
This is equivalent to derive the rotation curves assuming a con- 
stant CTgas ~ 20kms~-'^ for the halo region. 
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Figure 7. Rotational velocities of the model 6 = 0.1 and k = 1 
(solid line) and rotation curves in the HI by Fraternali et al. 
(2005) (points with error bars) at z = 3.9 kpc (top) and z = 5.2 
kpc (bottom). For each plot the long-dashed line represents the 
rotational velocity on the mid-plane. 



20kms~^ only at ^ 10 kpc (Fig. |4|. Ideally one should cor- 
rect the data points for this effect obtaining slightly lower 
velocities for R < 10 kpc and higher for R > 10 kpc but the 
effect is minimal. Overall, for what concerns the rotation 
curves of NGC 891 we can say that the high-cr^^ stationary 
model predicts a gradient that is quite close to (although 
slightly lower than) the observed one. 

The stationary models presented here do not aim to ex- 
plain the origin of the extra-planar gas. Instead their goal 
is to find a stationary description for the extra-planar gas 
phenomenon whatever physical mechanism is producing it, 
assuming that it can be considered roughly stationary. This 
assumption should be reasonably valid if the extra-planar 
gas is produced for instance by a galactic fountain and the 
star formation rate of the galaxy does not vary considerably 
with time as it seems to be the case for galaxies like the 
Milky Way (e.g. Twarog 1980). The investigation presented 
by Fraternali & Binney (2006) shows that the rotational ve- 
locities predicted by a galactic fountain model are system- 
atically higher than the one derived from the data and also 
higher than those predicted by the stationary high-cr^^ model 
shown in Fig. [71 This latter result is not surprising given 
that in the fountain model the gas clouds move outwards 
for a large fraction of their trajectories and thus tend to be 
at V(^ > Uc for most of their time (uc is the circular speed 
at that position). On the contrary, in the high-cr^ model 
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the gas has basically ^ Uc (Fig. |4]) . Fraternali & Binney 
(2008) revised the fountain model by including interactions 
of the fountain clouds with the surrounding ambient medium 
(see also Barnabe et al. 2006 for a similar suggestion). These 
interactions cause the fountain clouds to lose part of their 
angular momentum and to decrease their rotation veloci- 
ties. Such a model is to date the best description for the 
dynamics of the extra-planar gas. In addition to the rota- 
tional velocities, this latter model reproduces the data cube 
in detail including the channel maps close to the systemic 
velocity (see Fig. 2 of Fraternali & Binney 2008). Kaufmann 
et al. (2006) also reproduced the rotation curves of NGC 891 
fairly well with a model in which the extraplanar gas forms 
by condensation of the extended warm-hot corona, but they 
did not attempt to build pseudo-data cubes to compare 
directly with the observations. However, the thermal insta- 
bility on which such a model is based is unlikely to occur in 
galactic coronae, because of the combined effect of buoyancy 
and thermal conduction (Binney, Nipoti & Fraternali 2009). 

Whatever the origin of the extra-planar gas, it is im- 
portant to investigate whether or not it is in a stationary 
state. If so, it should be possible to find stationary models 
that reproduce these data in detail. This paper has been the 
second attempt of this sort after Barnabe et al. (2006). In 
both these attempts non-circular motions of the gas clouds 
have been neglected for simplicity, but in the kinematics of 
the extra-planar gas this kind of motions are clearly present 
(e.g. Boomsma et al. 2008; Fraternali et al. 2001) and they 
could affect substantially the observed kinematics of the 
halo. More likely, the cause of the discrepancy between our 
treatment and the observations lies on the assumptions of a 
Jeans-based model and in particular on the ballistic nature 
of the cloud motion. Not only the drag between the HI halo 
clouds and the corona may have an important dynamical 
effect but also at the turbulent cloud/corona interface there 
may be a considerable exchange of mass. Modelling this in- 
terface is a very challenging numerical problem (e.g. Vieser 
& Hensler 2007 a,b) and a secure estimate of the rate of mass 
exchange is not yet available. This issue will be studied in a 
forthcoming paper (Marinacci et al., in prep.). 



6 SUMMARY AND CONCLUSIONS 

In this paper, motivated by the results of Barnabe et al. 
(2006), we investigated the possibility that the extra-planar 
gas in spiral galaxies can be modelled as a "gas" composed 
by cold HI clouds that follows the stationary Jeans equa- 
tions. In doing this we are assuming that the pressure- 
bound clouds are moving almost ballistically in the halo. 
This assumption is valid in the limit of massive clouds hav- 
ing negligible rates of mass exchange with the corona. 

In this alternative interpretation of fluid baroclinic 
models the thermal pressure is replaced by an isotropic ve- 
locity dispersion tensor, so that the problem of high tem- 
perature of gaseous homogeneous models is eliminated. We 
have also extended the discussion to simple phenomenolog- 
ical models with anisotropic velocity dispersion tensors. We 
constructed both isotropic and anisotropic models in the 
well- constrained gravitational field of the spiral galaxy NGC 
891, and compared their predictions to the observed kine- 
matics of the extra-planar gas in that galaxy. For each model 



we built a pseudo-data cube with the same resolution and 
total flux as the observations. The main results of our anal- 
ysis can be summarized as follows: 

(i) The adopted functional form of the cloud density dis- 
tribution, taken by the observations, leads to physically ac- 
ceptable solutions in all the models investigated. The cloud 
density distribution is centrally depressed, and, in order to 
match the vertical extension of the HI halo of NGC 891, 
it has higher scale-height than that used by Barnabe et al. 
(2006). 

(ii) All the models computed show a negative vertical gra- 
dient in the rotational velocity, the distinctive feature of the 
kinematics of the extra-planar gas. 

(iii) The support against the vertical gravitational field 
of the galaxy requires a cr^ ^ 50 — lOOkms"^ and therefore 
the Une-of-sight velocity dispersion of the isotropic model is 
a factor ~ 3 — 4 higher than that observed. 

(iv) With the introduction of the anisotropy it is pos- 
sible to restrict the line-of-sight velocity dispersion to the 
observed values, but the predicted vertical gradient in the 
rotational velocity is somewhat too shallow and other fea- 
tures of the data cube are not fully reproduced. 

We conclude that the dynamics of extra-planar gas in a 
galaxy like NGC 891 is not fully described by any of the sta- 
tionary models considered here. However a model with an 
anisotropic velocity dispersion tensor, which mimics a galac- 
tic fountain is the preferable among all. The fact that none 
of the stationary Jeans models analysed here can reproduce 
all the features of the observed (extra-planar) gas kinematics 
might suggest that the cloud motion is not purely ballistic, 
and that the interaction between the clouds and the coronal 
gas, perhaps in the form of mass exchange, plays an impor- 
tant dynamical role. 
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